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1. INTRODUCTION 

During the past few decades, the fractional calculus gain much consideration because of its 
applicability in a lot of science and engineering fields like: electrical networks, phenomena in the fluid flow 
areas, probabilistic and statistics based decisions, chemical-physics, electrochemistry, and signal processing. 
That generally modelled by linear/nonlinear fractional differential equations [1]—[3]. To solve equations of 
different categories like linear or nonlinear, ordinary differential or partial differential equations, integer or 
fractional a number of methods have been used, for example adomian’s decomposition method [4]-[6], 
homotopy perturbation method [7], he’s variational iteration method [8], homotopy analysis method [9]. 

Among a multitude of available techniques to tackle nonlinear equations, homotopy analysis method 
(HAM) has gained tremendous popularity. Initially, Liao proposed in his doctoral thesis [10]-[13] the main 
method to develop the basic concepts of homotopy analysis in topical geometry to propose a general method 
of analysis for solving nonlinear problems. The homotopy analysis method does not depend on a small 
parameter, unlike other analytical techniques such as the method of decomposition, the homotopy 
perturbation method (HPM) as these methods are particular cases of the homotopy analysis method and 
which gives us a sufficient way to control and modify the convergence of solution chain by controlling the 
value of error by means of the convergence control parameter h. In this chapter, we have used homotopy 
analysis method with Padé approximations to solve linear partial differential equations of fractional order 
[14]-[16], the general formula of which is: 


ðu 


ate 


+ AC) + BO) = + Cu = h(x, t) (1) 
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where 
(x,t) € [0,1] x [0,T], n—-1<a< nneN (2) 


and initial conditions 


ky 
a(x, 0) = f(x) k=0,1,...,m-1 (3) 


2. DEFINITION 

Let f(x),x > 0 be a real function which is in the space CM, MeR where the real number k > M, 
for f(x) = x kg(x) where g(x) € [0,0) which is in the space CmM if and only if f(x)” E€ Cum EN 
[17]. 


2.1. Riemann-Liouville fractional integer 
The Riemann-Liouville fractional integral operator of order of a function f(x) € Cp, u 2 —1 
defined as in (4) and (5) [18]: 


IFO =F @-O YF Oadt,x> 0 (4) 

PF @) = f) (5) 
for (x) E€ Cuyu =—1,a,6 =0,y = —1 properties of the operator Ja 

LPG) = JE JEF) (6) 

jea = ee o 


2.2. Devrative Caputo fractional 
The fractional derivative of (x) [19] in the Caputo sense is defined as (8) [20]-[22]: 


DEF) = JD" FA) = ah @ ye F™ Wat (8) 


for (x) E (uw =-1,a,B = 0,y = —1,n—1 < æ < n,n E N properties of the operator Da 


DE DE FG) = Dy! f) = DE DEF GO (9) 
DŽ xY = we xt x > 0 (10) 
3. METHOD 


Let’s get the following nonlinear differential equation that is formed [23]-[25]: where N (nonlinear 
operator) (unknown function) (the independent variable). Let yọ be a first guess of the analytical solution y, 
h + O an assistant parameter, H(t) # OV t E R an assistant function and L an assistant linear operator with 
the property that L[y(t)] = 0 when y(t) =0. Then using q € [0,1] as parameter for embedding. We 
construct a homotopy that's called zero-order deformation: 


N[y@)]=0 fort =0 (11) 


Having considerable selection of the initial guess yọ should be stressed, the assistant linear operator. Here is 
L, and nonzero assistant function H(t). When q=0, can see the deformation equation zero-order (12) becomes 


(1 — q) LIO(t,q) — yo(t)] = qhH@)NI[O(t,@ )] (12) 


and when q=1, “since h 0 and H(t) /= 0”, the zero-order deformation (12) is tantamount to 
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O(t,0) = yo(t) (13) 
and when q = 1, since h + Oand H(t) + 0, the zero-order deformation (12) is equivalent to 

O(t.1) = y(t) (14) 
Thus, based on (13) and (14), the embedding-parameter q rises from 0 to 1, O(t,q) Different continually 


from y9(t) to the same solution y(t). This type of continous variation is deformation in homotopy. 
Also, Taylor’s theorem provides, @(t, q) might deal with a power series of q in the same expand: 


O(t,q) = Yo(t) + Lines Vm) q” (15) 
where 
a" Ot, 
Y(t) = SO Igo (16) 


If yo, L, h and the power series (9) of (t, q) converge at q = 1, then we have the solution series under those 
assumptions: 


y(t) = OCt 1) = yo(t) + Linz Ym) (17) 
for brevity, define the vector: 

In) = (¥0(t), 10), 20), nO} (18) 
According to (17), the governed-equation of y,,(t), also derived using the zero-order as in the (12) taht 


accour by the zero order deformation differentiation in (12) m -times with respective to q later to be divided 
by m! and later to set q = 0, to get m th -order deformation: 


L[Ym(t) -XmYm-1(t)] = hH (t) Rm-1 Omi) (19) 
where: 
> a™-1n[O( tq) 
Ra Fm) = a Nano (20) 
and 
Xm = {0, mad, 1, m>1 (21) 


4. PROPOSED ALGORITHM 

The algorithm was rewritten the homotopy analysis method and then linked to the pseudoPadé 
approximations. To find the initial evalution u0(x, t), we integrate (1) by Riemann-Liouville integration 
method with respect to the variable t and substituting the initial conditions to obtain 


u(x,t) = a w(x, 0%) + JE (AG) + Bo) HE + Cu) +h) (22) 


By omitting the term JaA(x)du +B(x) 02u+C(x)u in the right side in (22), we obtain the initial evolution 
u(x, t) by the formula 


k 
ug (x, t) = Eriu (x, 0+) — + J7h(x, t) (23) 


By using the deformation (19) of order m with assuming L = Da, H(t) = 1 where n E N and n-1 <a <n, 
we obtain the following sequential formula: 


UD“ Um (x,t) = Xm D“ Um_1 (x,t) + h Rm-1 (eee t)) (24) 
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aulro) a ABE) 


Rm- (tna) = 2 ( DOG, tq) + AG) + Bx) 


C(x)u(O(,t,q)) — AG, t) ) lq=0 (25) 


To find approximate iterations u,(x,t), u2(x,t), u3(x,t), we apply Ja Riemann-Liouville inte- gration on 
The both side of (24) to obtain 


: j 
Ua (Xs t) = Xmum-1(x, t) — Xm ‘ary Um-1” (x, 0+) T + hJE| Rm-1 (um-1(x, t))] (26) 


By substituting values of m=1, 2, 3, ... in (26), we obtain the approximated solutions u,(x,t), u2(x, t), 
u3(x, t). By substituting the approximated iterations u4 (x, t), uz (x, t), u3 (x, t),with ug(x, t) in the solution 
series Øm (x, t), we get: 


C6502 Zt us (x,t) (27) 


Then, substituting the values of h and a, we obtain an approximated solution for initial values problems (27) 
given by the form: 


Øm (x, t) = Ef-0 ath (28) 


where ao, a1,..., ak are constants, Bo, Bi,..., Bx are different positive powers, m is the number of iterations and 
k is specified according to the given example. 
Now, we connect the series (28) by Padé approximation using the following assumption: 


n-1 


ð u 
u(x,t) = >. u(x, Da +JE (ac 5 + BO) =a 


k=0 


t? = (29) 


2 


a cou) + J*h(x, t) 


where œ is fraction number, after substituting in the approximated solution (28), we obtain a series that has 
the form 


ea = Ek- Me ha (30) 


where Vo, Vi, ..., Vk are positive powers. By using Padé approximation of order M and N where M, N E N U 
{0}, we can convert the sequence (30) to a fractional series: 


Pn(z) 
Ry (2) = Gace) G1) 


Finally, by substituting the transformation z = tin the series (31), we obtain the Padé series with fractional 
powers. 


5. NUMERICAL EXAMPLES 
Suppose we have the following one dimensional heat equation that has a fractional order [4], [19], 


[20]. 
Du =} Uxx:0 <a< 1fort>0 (32) 
With condition (x, 0)=sin nx, The exact solution for the (32) is 


(-t%* 
T(aK+1) 


u(x,t) = sin nx VR-o (33) 


The solution: depending on the algorithm mentioned in the fifth paragraph, we get approximations: 
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U(x, t) = sin mx 


2xht* Fi 2ht2% $ 2ht* 
T(a+1) T(2a +1) Ta@+1) 


u(x,t) = 


ht*sinnx  h*t?“cosnx 


ule = (Veh) aay? Teak) 


We repeat this process to obtain u,4(x,t). By adding approximated iterations uo(x, t), u(x,t), uz2(x, t), 
U3 (x, t), u4(x,t), we obtain an approximated solution Ø5(x, t) where: 


ht“sinnx 2h cosnx 


_ v4 ee 2 h2t 
Øs (x, t) = Viz Um (%, t) = sinmx + TERED) (3+5h+2h*) + “Ton (2+h)+ 
h?t?%sinnx n3t3%cosnx 
+(1 + 2h) “Tara. +(1+h) EEN (34) 


Now, by taking different values of a with the best value of h at each value of a in the solution series (34) and 

then relating them with Padé approximation: 

— When o=0.25 and h=-0.18 and assuming that t°?5 = z, then calculating [1/2] the Padé’s 
approximations, and remembering that z =t°5, we get: 


sin(3.1415x) — 1.1238t°25sin(3.1415x) 
1 —.5.1442e71t9-25 4 4.794005 


[1/2] = (35) 
— When a=0.5 and h=0.25 and assuming that t°> = z in the solution series@ 5(x, t), then calculating [1/2] 


the Padé approximations, and remembering that z=t°>, we get: 


sin(3.1415x) — 1.4094t955'"(3,1415x) 
1 — 6.3806,°° + 2.2221t 


[1/2] = (36) 
— When o=0.75 and h=0.22 and Assuming that t°75 = z in the solution series 05(x, t),then calculating 


[1/2] the Padé’s approximations, and remembering that z = t°-75, we get: 


sin(3.1415x) — 1.2292t®-75 sin(3.1415x) 
1-—5.4392e71,°°7° + 1.3873t15 


[1/2] = (37) 


Suppose we have a one dimensional linear no homogenous Burgers equation of fractional order 


[26]. 

Diu + Uy — Ur = + 2x 2 (38) 
when: 

0<a<1,t0,x E€ R With condition u(x, 0) = x? (39) 


the exact solution for the problem (38) is 
u(x,t) = x? + t? (40) 


The solution: depending on the algorithm mentioned in the fifth paragraph, we get approximations: 


a 


Ug (x, t) = x? + t? tiero = 2) 
jz 2xht® $ 2ht2% 2ht® 
m0 =F aa treat Tat 
ye 2xht®” 2ht* 5 2ht2* j 
u(x,t) = (1 + 2h +h?) (1+h+h?) + (1 + 4h + 3h?) 


T(a+1) T(a@+1) T(2a +1) 
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We repeat this process to obtain u,(x,t). By adding approximated iterations ug(x, t), u,(x,t), u(x,t) , 
U3 (x,t), u,(x, t), we obtain an approximated solution @5(x, t) where: 


D =x? +t? + —_ (2x- gy eat. 
Os (x,t)= Ef -0 Un (x t)=x? +t Faan (2x-2 )+(3+5h+2h lean 
2ht® 2 2ht?* 2113 
ei (3+3h+2h ETOT (3+10h+8h?+h3)+ (41) 


Now, by taking different values of a with the best value of h at each value of a in the solution series (41) and 

then relating them with Padé approximation: 

— When o=0.25 and h=—0.7 and assuming that t°?° = z in the solution series Ø5(x, t), then calculating 
the PA [5/6] the Padé’s approximations, and remembering that z = t°°, we get: 


ae x2 + e-?9(1.7872e22x —1.7872e27)t0.25 —1.7061e — 1t2(42) 
6 (1 + (34.4190x2 — 1.2864e — 1x + 6.4321e — 2)t + (6.1401e — 1 — 6.1401e — 1)t1.25 + 5.8612t15) 


e-2963.4419e3%x* 1.2864e78e3+6.4323e77x2t+e~7941.2291e 29x? 1.2314€79x73.4489H, 
(:14+W(3h4e.4n190x? —1.2864e1 x+6.4321e~2 )t+(6.1401e-1-6.1401e71 )t1-25 +5.8612t15) 


(42) 


— When a=0.5 and h=0.85 and assuming that t°° = z in the solution series 5(x, t),@ then calculating [1/2] 
the Padé’s approximations, and remembering that z = t°”, we get: 


3.5861 e°4x444.0789e°0x3+2.0394e90x2 
E ETET] 
(+ 4.0968e°1x—4.0968e°1 ) 0.54. 8.2297e52t ) 
3.5861e°4x2—4.0789e29+42.0394e90 3.5861e°4x2—4,0789e20+42.0394e90 


PA | = (43) 


— When o=0.75 and h=0.6 and Assuming that t°75 = z in the solution series Ø5(x, t), then calculating 
[2/3] the Padé’s approximations, and remembering that z = t°”°, we get: 


x2 


3.2221e— 56(-1.7289e94x+1.7289e94)t0.75 
1+ my. 


PA[2/3] = (44) 


6. RESULTS 

The results of this paper show that the proposed homotopy method analysis solve partial differential 
equations efficiently. The results have been verified with other works to show that the proposed method has a 
better mean square error and exact solution. The results compared with both HAM and HAM-PA when 
t=0.01 and 0.1 as shown in Tables | to 5. 


Table 1. The obtained results which compared to the Table 2. The obtained results which compared to the 


exact solution for instance example 5.1, at t=0.01 exact solution for instance example 5.2, at t=0.01 
x Exact Solution HAM HAM-PA x Exact Solution HAM HAM-PA 
0 0 0 0 0 e* —2.2612e7?  —5.6468e-3 
0.1 2.5779e —1 2.5518e — 1 2.5639e — 1 0.1 1.01e7? —1.2047e7? 4.8589e7? 
0.2 4.9035e — 1 4.8538e — 1 4.8769e — 1 0.2 4.01e7? 1.851e7? 3.5250e7? 
0.3 6.7491e — 1 6.6807e — 1 6.7124e — 1 0.3 9.01e7? 6.9082e7? 8.5535e7? 
0.4 7.9340e — 1 7.8537e —1 7.8910e — 1 0.4 1.601e7+ 1.3964e7+ 1.5572e77 
0.5 8.3423e — 1 8.2578e — 1 8.2970e — 1 0.5 2.501e77 2.3021e77 2.4582e77 
0.6 7.9340e — 1 7.8537e— 1 7.8910e — 1 0.6 3.601e77 3.4077e77 3.5586e7* 
0.7 6.7491e — 1 6.6807e — 1 6.7124e — 1 0.7 4.901e7t 4.7134e71 4.8585e7t 
0.8 4.9035e — 1 4.8538e — 1 4.8769e — 1 0.8 6.401e77 6.2190e7+ 6.3579e7* 
0.9 2.5779e —1 2.5518e — 1 2.5639e — 1 0.9 8.101e71 7.9247e7! 8.0572e71 
1.0 —4.1471e — 31  —41051le—1 —4.1246e — 1 1.0 1.0001 9.8303e77 9.9564e71 
Table 3. Comparison of the mean square error for Table 4. Comparison of the mean square error for 
each of the HAM and HAM-PA [5%] each of the HAM and HAM-PA [1/2] 
Algorithm Name Error Algorithm Name Error 
MSE- HAM 3.9868e — 4 MSE-HAM 8.3458e — 8 
MSE - HAM-PA 2.1430e — 5 MSE - HAM-PA 7.7009e — 8 
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7. 


Table 5. Comparison of the mean square error for 
each of the HAM and HAM-PA [2/3] 


Algorithm Name Error 
MSE-HAM 1.5467e — 6 
MSE - HAM-PA 8.7674e — 7 


CONCLUSION 
This study, the efficiency of the HAM with Padé approximations to solve partial fractional linear 


questions has been proven. The results showed the possibility of Padé approximations to improve the results 
of the HAM. Padé approximation only solves sequences with natural bases. Therefore, we used the 
hypothesis to=x where œ is the number of fractional basis to normal basis. 
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